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Abstract 

The  effect  of  self-heating  and  cooling  by  natural  convection  on  a  sustainable  temperature  of  PEM  fuel  cell  stacks  was  studied.  Overall  mass 
and  heat  balance  equations  are  combined  to  predict  self-heated  temperatures  at  various  operating  conditions.  Analyses  show  that  the  effect 
of  a  heat  loss  coefficient  is  more  important  than  other  variables  such  as  air  flow  rate  and  surrounding  temperature.  The  stack  design  variables 
such  as  active  cell  area  and  number  of  cells  also  have  significant  influence  on  self-controlled  temperature.  A  lower  Ohmic  resistance  of  cells 
is  expected  to  allow  a  wider  range  of  current  load  applications.  The  proposed  model  can  also  be  used  to  evaluate  heat  loss  coefficient  from 
measured  stack  performance  and  temperature  data.  Experiments  performed  on  a  seven-cell  stack  of  50  cm2  active  area  were  used  to  provide 
data  for  the  validation  of  the  model. 

©  2005  Elsevier  B  .V.  All  rights  reserved. 
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1.  Introduction 

Proton-exchange-membrane  fuel  cells  (PEMFC)  have  a 
narrow  range  of  operating  temperatures  due  to  the  need  of 
water-containing  electrolyte  and  the  limit  of  low  tempera¬ 
ture  for  proton  conductivity.  Temperature  control  or  thermal 
management  is  a  practical  issue  in  general  for  the  design  of 
fuel  cell  systems  [1].  Unlike  the  high  temperature  fuel  cells,  a 
PEMFC  works  at  a  relatively  low  temperature  and  it  does  not 
have  to  rely  heavily  on  the  use  of  external  thermal  manage¬ 
ment  systems  for  heating  or  cooling.  At  a  temperature  below 
100  °C,  a  self-controlled  temperature  is  conceivable  by  bal¬ 
ancing  heat  released  from  the  membrane-electrode-assembly 
(MEA)  and  heat  loss  through  end-plates  to  surrounding  at¬ 
mosphere.  Self-heating  makes  the  system  simpler,  just  as 
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self-humidification  does.  While  self-humidification  has  been 
mostly  a  material  design  issue,  self-heating  and  cooling  is  an 
engineering  design  issue. 

Both  the  heat  release  and  heat  removal  rates  are  depen¬ 
dent  on  many  variables.  The  objectives  of  this  study  are:  (1) 
to  identify  these  variables  and  introduce  a  model  that  ac¬ 
counts  for  those  aspects  of  self-heating  and  cooling,  (2)  to 
demonstrate  how  such  a  model  can  be  applied  to  the  design 
of  stacks.  An  overall  mass  and  energy  balance  model  is  devel¬ 
oped  in  this  work.  Although  we  focused  on  using  the  model 
to  predict  temperature  of  PEM  fuel  cells  in  the  absence  of 
external  heating  and  cooling  devices  in  this  work,  the  model 
developed  here  can  be  applied  to  any  kind  of  fuel  cells  in 
general. 

2.  Overall  heat  and  mass  balance 

We  propose  a  model  based  on  a  combination  of  the  steady- 
state  overall  heat  and  mass  balance  equations.  With  hydrogen 
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coming  in  through  anode  and  air  through  cathode,  both  of 
humidified,  the  overall  mass  balances  for  each  species  are 
expressed  in  mass  or  molar  flow  rates  (nj)  and  the  reaction 


rate  (rn2o)  by  Faraday’s  law. 

^H2,in(A)  -  m20  =  ^H2,out(A)  (1) 

^H20,in(A)  =  ^H20,out(A)  (2) 

^02,in(C)  -  0-5rH2O  =  ^02,out(C)  (3) 

^N2,in(C)  =  ^N2,out(C)  (4) 

^H20,in(C)  +  m20  =  ^H20,out(C)  (5) 


Water  is  assumed  to  be  present  in  the  saturated  vapor  phase 
only,  for  simplification  of  the  model.  A  saturated  water  va¬ 
por  pressure,  Psw,  can  be  calculated  from  the  saturation  tem¬ 
perature,  Tsw,  using  the  thermodynamic  equation  (Antoine 
equation): 


1ft  P sw 


18.3036  -  3816.44 
Fsw- 46.13 


where  the  unites  of  P sw  is  in  mmHg  and  Tsw  is  in  Kelvin. 

Humidified  gas  feed  rates  are  calculated  with  the  added 
saturated  water  vapor  as: 


Fsw 

P -  =  ^H2o 

*  system 

^H2Q(A)  _  yH2Q(A) 

^H2(A)  1  -  yH20(A) 

^h2o(C)  =  yH2o(C) 

^o2(C)  +  ^n2(o  l  -  yu2o(C) 


(7) 

(8) 
(9) 


The  heat  balance  is  obtained  from  the  enthalpies  of  the 
anode  and  cathode  gas  streams  at  inlet  and  outlet  with  the  ref¬ 
erence  temperature  (To,  usually  25  °C).  Self-heating  is  calcu¬ 
lated  from  heat  generated  during  the  fuel  cell  reaction  (Gcell)- 
Self-cooling  is  calculated  from  heat  loss  to  an  atmosphere 
(Gloss)  with  a  surrounding  temperature  (Tsurr)  through  a  to¬ 
tal  external  surface  area  of  exposure  (Aext)-  External  heating 
or  cooling  is  not  considered  in  this  model.  The  heat  balance 
equations  are: 


input  enthalpy  =  Win(A)C/?,in(A)(7in(A)  -  To) 

+  ^in(C)6,p,in(C)(7in(C)  “  7o)  (10) 


output  enthalpy  =  ft0ut(A)Cp,out(A)(7but(A)  -  To) 

+  ^out(C)C/?,out(C)(7but(C)  —  To)  (11) 

heat  generated  =  Gcell  =  (-ATT)H20rH2o-^ceU  ^2  yce11 

(12) 

heat  loss  =  Gloss  —  h loss  Aext (Tcell  Fsurr)  (13) 

The  heat  capacity  of  gas  mixture  is  easily  obtained  from 
those  of  pure  gases,  and  the  enthalpy  change  of  fuel  cell 


reaction  is  shown  in  thermodynamic  tables  as  a  function  of 
temperature  [2]. 

6^?,  mix  (y  j  i  T)  =  J2yjcpj(T)  and  AH  =  f(AH0,T ) 

(14) 

There  is  a  temperature  distribution  inside  a  fuel  cell  stack. 
In  the  present  simplified  model,  the  variation  of  temperature  is 
neglected  and  a  single  average  temperature  (Tcen)  is  assumed, 
with  a  further  assumption  that  outlet  gas  temperature  is  equal 
to  this  average  cell  temperature. 

Tcell  —  W)  ~  7but(C)  (15) 

Combining  Eq.  (10)  through  (15)  and  the  assumptions 
made,  we  obtain  the  following  simple  equation  for  tempera¬ 
ture  predictions  of  self-heated  fuel  cells: 

^in(A)6’p,in(A)(7in(A)  To) 

T  ^in(C)6"p,in(C)(7in(C)  —  To)  T  Gcell  —  Gloss 

—  (^out(A)Cp?0ut(A)  T  ^out(C)6’p,out(C))(Tcell  —  To)  (16) 

The  overall  balance  model  and  parameters  are  schemat¬ 
ically  shown  in  Fig.  1.  Sa  and  Sc  are  reactant  stoichiome¬ 
try  (defined  here  is  a  ratio  of  the  reactant  gas  flow  rate  di¬ 
vided  by  the  amount  required  for  100%  conversion  to  gen¬ 
erate  electrical  current  applied)  of  anode  and  cathode  gases 
for  a  given  current  load.  Thermodynamic  properties  (heat 
capacity  and  enthalpy  change  of  the  reaction)  are  depen¬ 
dent  on  temperature.  The  enthalpy  change  of  fuel  cell  re¬ 
action  (H2  +  I/2O2  — >  H2O)  is  represented  by  a  linear  cor¬ 
relation  presented  in  Table  1.  There  is  only  a  0.2%  varia¬ 
tion  of  the  enthalpy  change  for  a  50  °C  temperature  change. 
Therefore,  we  assume  a  constant  enthalpy  change  over  the 
temperature  range  of  interest.  Such  a  negligible  depen¬ 
dence  on  temperature  is  also  observed  in  heat  capacities. 
Table  1  shows  linear  correlations  of  enthalpy  and  heat  ca¬ 
pacities  and  their  constant  values  at  the  PEM  FC  temperature 
range. 

With  these  further  simplifications,  the  input  variables  in 
the  overall  balance  equation  are  all  constants  except  the  cell 
voltage  (Vceii)-  Cell  or  stack  voltage  is  dependent  on  current 
load,  humidity,  temperature,  pressure,  water  flooding,  and 
so  on.  Fuel  cell  models  are  often  very  comprehensive  and 
account  for  all  these  effects.  As  a  result  of  this  comprehen¬ 
siveness,  the  models  proposed  are  very  diverse  as  reported 


Table  1 


Linear  correlation  (A  +  BT)  and  average  values  of  thermodynamic  properties 


Parameter 

A 

B 

Avg.  (300-400  K) 

-Atf(H20) 

238870.3 

9.92008 

242346  J  mol- 1 

Cp( H2) 

28.209 

0.0022 

29.02  J  mol-1  K”1 

Cp(  o2) 

26.847 

0.00839 

29.75  J  mol-1  K-1 

Cp( N2) 

28.451 

0.00218 

29.19  Jmor1  K-1 

Cp(  H2o) 

31.171 

0.00801 

33.93  J  mor1  KT1 
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Fig.  1.  Overall  heat  balance  model  at  steady-state. 


in  the  literature  [3].  Because  so  many  variables  affect  cell 
performance  in  very  different  patterns  and  also  because  there 
is  no  single  equation  to  account  for  such  a  relationship,  to 
avoid  unnecessary  complication,  we  choose  the  well-known 
model  using  a  Tafel  slope  ( b )  and  Ohmic  resistance  (Rohm) 
to  calculate  the  cell  voltage: 

^cell  =  Eo  —  b  In  i  —  iRohm  (17) 

3.  Experimental  setup 

A  commercially  available  small  PEM  fuel  cell  stack  was 
used  to  provide  experimental  data  for  model  validation.  The 
stack  has  seven  MEAs  (Nafion-115,  lmgcm-2,  20wt.% 
Pt/C)  with  an  active  area  of  50  cm2  (EFC-50-ST,  Electrochem 
Inc.,  Woburn,  MA).  The  six  graphite  bipolar  plates  have  a 
dimension  of  4  in.  x  4  in.  x  1/4  in.,  and  two  graphite  end- 
plates  have  a  dimension  of  4  in.  x  4  in.  x  1  /  2  in.  with  the 
two  end  frames  of  5  in.  x  5  in.  x  1/2  in.  The  total  external 
surface  area  is  820  cm2.  We  tested  this  stack  in  a  test  system 
equipped  with  bubble  humidifiers  and  heated  gas  transfer 
lines  (FCTS-HTK  and  FCTS-HB,  Lynntech  Industries  Ltd., 
College  Station,  TX).  Several  tests  were  done  at  various  cur¬ 
rent  loads  and  humidifying  conditions,  as  shown  in  Table  2.  In 
each  test,  current  load  increased  from  zero  to  a  certain  value, 
and  the  temperature  rise  and  the  change  of  stack  voltage  were 
observed  until  they  reach  steady-state  values. 

4.  Results  and  discussion 

The  variables  in  the  model  schematics  shown  in  Fig.  1  are 
classified  into  three  categories:  controllable,  dependent,  and 


noise  variables.  Controllable  or  input  variables  are  those  we 
can  control  independently,  and  they  are  reactant  stoichiome¬ 
tries  (Sa  and  Sq),  saturated  water  temperatures  (Tsw),  heating 
line  temperatures  (7A),  etc.  Active  cell  area  (Acen),  number 
of  cells  for  stack  (Abells),  and  external  area  of  exposure  (Aext) 
are  also  controllable  at  the  design  stage,  though  not  during 
operation. 

To  examine  the  effects  of  these  variables  on  self- 
maintained  temperature  in  PEM  fuel  cells,  calculations  were 
done  for  the  seven-cell  stack  (described  in  Section  3)  at  a  bub¬ 
ble  humidifier  temperature  of  75  °C,  a  surrounding  tempera¬ 
ture  of  20  °C,  and  a  heat  loss  coefficient  of  10.0  W  m~2  K_1 . 
The  relationship  between  current  density  and  cell  perfor¬ 
mance  (voltage)  needs  empirical  parameters,  and  we  chose 
the  reference  data  for  the  empirical  parameters  related  to 
the  Tafel  equation  from  the  literature  [4]:  E$  =  0.942  V; 
b  =  0.061  V;  and  ^ohm  =  0.39  £2  cm2.  These  are  for  Nafion 
115  at  50  °C. 

4.1.  Effect  of  operating  condition  and  design  variables 

Fig.  2  shows  the  effect  of  one  of  the  controllable  vari¬ 
ables,  air  stoichiometry.  Self-heated  temperature  was  cal¬ 
culated  by  changing  air  stoichiometry  from  1.0  to  10.0 
while  other  variables  are  unchanged.  The  fuel  (hydrogen) 
stoichiometry  was  also  kept  constant  at  1.5.  In  fuel  cells, 
the  air  stream  to  cathode  is  frequently  used  for  the  pur¬ 
pose  of  cooling,  especially  in  large-scale  high-temperature 
fuel  cells  [5].  In  PEM  fuel  cells,  a  high  flow  rate  of  air  is 
also  useful  to  minimize  water  accumulation  in  the  cathode 
channel. 

Fig.  2  shows  that  the  high  flow  rate  of  air  (with  a  higher 
stoichiometry)  has  some  cooling  effect  at  a  relatively  high 
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Table  2 


Experimental  conditions  for  the  stack  tests  to  measure  self-heated  temperature 


No. 

Stoichiometry,  5 a AS’c 

Tsw  (°C),  H2/air 

1(A) 

Vstack  (V) 

Favg  (°C) 

hioss  (Wm-2  K->) 

1 

1.20/1.20 

0/0 

1.0 

1.988 

28.06 

13.64 

2 

1.20/1.20 

0/0 

2.0 

1.346 

29.58 

23.86 

3 

1.20/1.20 

0/0 

3.0 

0.711 

36.29 

20.54 

4 

2.82/2.96 

0/0 

4.0 

0.694 

43.65 

17.75 

5 

2.82/2.96 

50/50 

4.0 

1.09 

47.06 

15.08 

6 

2.82/2.96 

22/50 

4.0 

0.813 

46.98 

15.50 

7 

2.25/2.37 

22/22 

5.0 

0.735 

50.77 

16.57 

8 

2.25/2.37 

50/22 

5.0 

1.176 

51.69 

15.27 

9 

3.01/3.16 

70/60 

5.0 

2.324 

57.6 

11.24 

10 

3.76/4.73 

53/53 

2.5 

2.0 

43.0 

10.09 

current  load.  At  a  relatively  low  current  load  of  0.25  A  cm  , 
the  air  cooling  has  almost  no  influence  on  temperature.  At 
such  a  low  current  density,  the  consumption  of  fuel  is  very 
small  and  the  rate  of  heat  generation  is  too  low  to  heat  the  fuel 
cell  stack.  The  temperature  at  0.25  A  cm-2  is  even  slightly 
increased  as  air  stoichiometry  increases.  This  is  because  the 
gas  inlet  temperature  (75  °C)  is  higher  than  the  fuel  cell  tem¬ 
perature.  The  increased  air  stream  carries  more  heat  from 
the  humidified  line  to  the  fuel  cell  stack  than  that  produced 
from  the  stack  itself,  as  the  heat  generation  rate  is  very 
low.  If  the  humidifying  condition  changes  to  a  lower  tem¬ 
perature,  we  will  see  the  air  cooling  effect  by  an  increased 
air  stoichiometry,  though  not  significant  at  a  low  current 
load. 

Figs.  3  and  4  show  the  effects  of  stack  size  in  terms  of 
active  area  and  number  of  cells.  In  single-cells,  there  is  al¬ 
most  certainly  no  effect  of  self-heating  on  fuel  cell  temper¬ 
ature.  Single-cells  are  mostly  used  for  fundamental  research 
purposes  in  laboratories.  In  laboratory  tests,  cell  tempera¬ 
ture  usually  needs  to  be  controlled  at  a  desired  set  point. 
Therefore,  single  tests  need  external  heating,  even  at  a  high 
current  load.  Larger-scale  stacks,  on  the  other  hand,  have  a 
distinct  self-heating  effect  resulting  in  a  significant  tempera¬ 
ture  rise.  More  often,  they  require  cooling  than  heating  unless 
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Fig.  3.  Effect  of  active  area  on  self-heated  temperature  at  a  current  density 
of  0.25  A  cm-2. 


the  current  load  is  extremely  low.  As  we  compare  the  cases 
of  Figs.  3  and  4,  the  five-cell  stack  cannot  maintain  temper¬ 
ature  below  100  °C  at  a  current  density  of  0.5  A  cm-2.  This 
indicates  a  self-heating  stack  can  be  used  only  for  low  current 
applications. 

Fig.  5  shows  a  contour  plot  of  temperature  on  the  plane 
coordinate  of  active  area  and  number  of  cells.  This  plot  il- 
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Fig.  5.  Temperature  contours  (°C)  for  the  effect  of  stack  size  (active  area 
and  number  of  cells)  on  the  self-heating  PEM  fuel  cell  stack. 

lustrates  a  design  issue  of  stack  dimension.  Even  at  a  low 
current  application,  a  self-heating  stack  is  to  be  designed  to 
have  at  least  one  stack  dimension  limited:  either  the  num¬ 
ber  of  cells  should  be  small  or  the  active  area  should  be 
small. 

4.2.  Effect  of  fuel  cell  performance  variables 

Dependent  variables  are  usually  output  variables  such  as 
cell  voltage,  which  is  strongly  dependent  upon  materials 
(membrane,  catalysts,  and  gas  diffusion  layers),  operating 
conditions,  and  hardware  such  as  cell  frames.  A  decrease 
of  cell  voltage  means  a  lower  efficiency  in  fuel  conversion  to 
electrical  energy  which  then  means  more  thermal  energy  pro¬ 
duced.  The  effects  of  two  cell  performance  parameters  (Tafel 
slope  and  Ohmic  resistance)  were  studied.  In  Fig.  6,  the  Tafel 
slope  varied  from  0.0 1  to  0.09  and  the  self-heated  temperature 
was  calculated.  Fig.  7  shows  the  effect  of  Ohmic  resistance 
when  it  changed  from  0.05  to  3.0  £2  cm2.  The  Tafel  slope  is 
a  measure  of  activation  overpotential  which  is  dominant  at 
relatively  low  current  load.  Due  to  the  relatively  low  amount 


Fig.  6.  Self-heated  temperature  as  a  function  of  current  density  at  various 
Tafel  slopes. 


Current  density  (A/cm2) 

Fig.  7.  Self-heated  temperature  as  a  function  of  current  density  at  various 
Ohmic  resistances. 

of  heat  produced  at  a  low  current  load,  the  effect  of  Tafel 
slope  is  not  significant  on  self-heated  temperature,  as  shown 
in  Fig.  6.  The  effect  of  Ohmic  resistance  is  more  significant 
because  it  is  the  dominant  factor  for  cell  performance  at  an 
increased  current  load.  Obviously,  sustainable  fuel  cell  tem¬ 
perature  is  obtained  for  a  broader  range  of  current  load  when 
ME  A  materials  have  lower  Ohmic  resistances,  as  shown  in 
Fig.  7. 

4.3.  Effect  of  noise  variables 

Noise  variables  are  those  which  have  significant  influences 
on  the  performance  variables  but  it  is  difficult  to  control. 
This  is  a  term  usually  encountered  in  statistical  analysis  or 
quality  engineering.  In  self-heating/cooling  PEM  fuel  cells 
and  stacks,  heat  loss  to  an  atmosphere  is  generally  difficult 
to  predict,  due  to  a  relatively  large  variance  of  surround¬ 
ing  temperature  (Tsurr)  and  uncertainty  in  a  heat  loss  coef¬ 
ficient  (Moss).  For  this  reason,  they  are  considered  noise  vari¬ 
ables.  We  can  reasonably  presume  that  the  surrounding  at¬ 
mosphere  temperature  varies  from  place  to  place  in  the  range 
of  0-40  °C.  The  heat  loss  coefficient  is  approximately  in  the 
range  of  3-12  W  m-2  K-1  for  natural  convection  through  the 
surrounding  air  and  in  the  range  of  12-85  Wm_2K_1  for 
forced  convection  through  air  [6].  A  total  external  area  of 
cells/stacks  exposed  (Aext)  can  be  a  noise  factor,  too,  because 
it  is  often  difficult  to  define  this  area  in  different  shapes  and 
configuration  of  cell/stack  frames  and  end-plates. 

Fig.  8  shows  that  a  high  current  load  operation  requires 
forced  convection  to  maintain  temperature  below  100  °C. 
With  only  natural  convection  near  room  temperature,  stack 
operation  is  limited  up  to  0.5  A  cm-2  in  this  specific  example. 
Forced  convection  usually  needs  a  device  with  an  external 
supply  of  energy  such  as  a  motor-driven  fan.  For  practical 
applications,  the  energy  to  drive  such  an  external  cooling  de¬ 
vice  must  be  provided  from  the  fuel  cell  stack.  Consider  a 
case  when  a  fuel  cell  device  needs  to  run  at  a  current  density 
of  0.75  A  cm-2  for  a  given  stack  design  and  performance: 
if  we  use  the  result  of  Fig.  8  for  further  design,  self-heating 
will  raise  the  stack  temperature  to  near  150  °C  with  natu- 
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Fig.  8.  Effect  of  the  heat  loss  coefficient  on  self-heated  fuel  cell  temperature 
at  various  current  loads  with  constant  stoichiometry  ratios  of  fuel  (Sa  =  1.5) 
and  oxidant  (Sc  =  2.5)  gases. 

ral  convection  cooling  (assuming  h\oss  =  10  W  m-2  K_1).  To 
maintain  temperature  below  80  °C,  the  stack  needs  a  cool¬ 
ing  fan  that  can  increase  the  heat  loss  coefficient  to  about 
30  W  m-2  K_1 .  It  will  take  a  couple  of  experimental  tests  to 
evaluate  how  much  power  is  needed  to  achieve  such  a  heat 
loss  coefficient.  If  the  expected  power  requirement  is  a  small 
portion  of  the  power  produced  from  the  fuel  cell  stack,  the 
design  of  this  self-heated  stack  with  an  auxiliary  cooling  fan 
is  feasible  with  respect  to  overall  heat  balance  and  power 
consumption.  This  case  illustrates  how  the  proposed  model 
can  be  quickly  and  effectively  used  for  design  of  fuel  cell 
application  systems. 

Fig.  9  shows  the  effect  of  another  noise  parameter,  sur¬ 
rounding  temperature,  on  self-heating.  The  surrounding  tem¬ 
perature  in  the  range  of  0-40  °C  does  not  have  a  very  signif¬ 
icant  influence  compared  to  the  other  parameters. 

4.4.  Evaluation  of  heat  loss  coefficient 

Table  2  summarizes  ten  experimental  tests  for  measure¬ 
ment  of  self-heated  temperature  at  various  conditions.  This 


Tsurr  (°C) 


Fig.  9.  Effect  of  surrounding  temperature  on  self-heated  fuel  cell  temper¬ 
ature  at  various  current  loads  with  constant  stoichiometry  ratios  of  fuel 
(, Sa  =  1.5)  and  oxidant  (Sc  =  2.5)  gases. 


seven-cell  PEM  FC  stack  showed  a  relatively  poor  over¬ 
all  performance,  and  it  was  difficult  to  apply  a  load  cur¬ 
rent  higher  than  5  A  (0.1  A  cm-2)  due  to  a  low  stack  volt¬ 
age.  Stack  voltage  was  unstable  during  the  constant  current 
load,  and  it  was  more  unstable  in  the  test  with  dry  fuel  and 
air  gases  (Tsw  =  0  °C).  Temperatures  were  measured  at  three 
different  points  (two  end-plates  and  one  bipolar  plate),  and 
they  were  averaged.  There  was,  in  fact,  a  small  degree  of 
temperature  variation.  The  inlet  end-plate  temperature  was 
lower  than  the  outlet  end-plate  temperature,  and  the  cen¬ 
ter  bipolar  plate  temperature  was  the  highest.  In  the  over¬ 
all  balance  model,  however,  an  average  temperature  is  suf¬ 
ficient  to  evaluate  a  heat  loss  coefficient  from  Eqs.  (10)  to 
(16).  The  resulted  heat  loss  coefficient  ranged  from  10.09  to 
23.86  Wm"2  K_1. 

The  results  in  Table  2  show  that  the  heat  loss  coeffi¬ 
cient  becomes  high  with  dry  gases,  and  also  it  becomes  low 
when  self-heated  temperature  increases.  Overall,  however, 
the  heat  loss  coefficient  does  not  make  a  reasonable  corre¬ 
lation  with  any  variables,  so  we  believe  this  heat  loss  coef¬ 
ficient  is  intrinsic  to  the  system.  It  will  vary  from  one  stack 
to  another,  depending  on  stack  size  and  materials,  but  it  is 
an  intrinsic  value  of  a  specific  stack;  therefore,  particular 
attention  is  required  to  the  heat  loss  coefficient  in  the  pro¬ 
cess  of  fuel  cell  stack  design  in  order  for  self-heating  to  be 
effecting. 

5.  Conclusion 

The  proposed  overall  balance  model  allows  us  to  study 
the  effects  of  various  parameters  on  self-heated  fuel  cell  tem¬ 
perature.  The  model  is  straightforward,  and  can  provide  an 
average  temperature  that  is  sufficient  for  the  purpose  of  a  ba¬ 
sic  design  of  specific  fuel  cell  application  systems  without 
external  heating.  From  our  case  studies,  we  can  draw  several 
findings  to  conclude:  (1)  self-heating  fuel  cells  can  be  effec¬ 
tive  for  low  current  loads  in  general;  (2)  the  design  of  stack 
size  is  important  and  at  least  one  of  the  stack  dimensions 
(area  or  number  of  cells)  should  be  small  for  self-heating  to 
be  effective;  (3)  the  range  of  allowable  current  load  increases 
with  lower  Ohmic  resistance  materials;  (4)  high  current  load 
applications  of  self-heating  stack  can  be  made  possible  by 
increasing  the  heat  loss  coefficient  to  the  range  of  forced 
convection;  (5)  heat  loss  coefficient  has  a  strong  effect  on 
self-heating  effectiveness.  We  also  found  from  experimental 
measurement  that  the  heat  loss  coefficient  by  natural  con¬ 
vection  is  mostly  in  the  range  of  10-20  Wm-2  K_1,  which 
agrees  with  the  typical  range  of  heat  transfer  coefficient  in 
atmosphere  known  from  the  literature. 
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